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Coarse-graining the constitutive equation 
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A recent theory for stress transmission in isostatic granu- 
lar and cellular systems predicts a constitutive equation that 
couples the stress field to the local microstructure [1]. The 
theory could not be applied to macroscopic systems because 
the constitutive equation becomes trivial upon straightfor- 
ward coarse-graining. This problem is resolved here for ar- 
bitrary planar structures. The solution is based on the ob- 
servation that staggered order makes it possible to couple 
the stress to a reduced geometric tensor that can be coarse- 
grained. The method proposed here makes it possible to 
apply this idea to realistic systems whose staggered order is 
generally 'frustrated'. This is achieved by a renormalization 
procedure which removes the frustration and enables the use 
of the upscalable reduced tensor. As an example we calculate 
the stress due to a defect in a periodic solid foam. 

Particulate and cellular materials play a major role in 
everyday life and are relevant to a wide range of tech- 
nological applications. Yet, the fundamentals of stress 
transmission in these systems are not fully understood. 
Recent work has shown that stress transmission in both 
cellular and isostatic granular systems differs significantly 
from conventional solids [2]- [5]. All solids in mechanical 
equilibrium satisfy force and torque balance 



V • CT i^ext = 



(1) 



where a is the stress tensor, cr^ its transpose and Ff.xt is 
an external loading. To the balance equations one must 
add a set of constitutive equations to uniquely determine 
the stress field. In elasticity theory these equations in- 
volve stress-strain relations but this information is redun- 
dant for granular [1]- [4] and cellular systems [5]. This is 
because these systems are isostatic and the intergranular 
and intercellular forces can be in principle determined 
from statics alone. Note that this is true regardless of 
material compliance. It has been shown that in 2D the 
conventional constitutive relation can be replaced by an 
equation of the form 



(2) 



with Q a tensor that depends only on the local geometry. 
Several empirical [3] and mean-field [4] proposals have 
been put forward for the macroscopic form of Q and a 
recent theory gave the specific dependence of this tensor 
on microstructural details [1]. In a further development 
an exact mapping was found between trivalent cellular 
solids and isostatic granular systems which enabled to 
extend the theory to solid foams [5]. For brevity, the 




FIG. 1. The vectors fcv and Rev shared between vertex v 
and cell c. rev connects two neighboring wall midpoints and 
is one edge in a clockwise-directed triangle around vertex v. 
The vector Rev points from the center of triangle v to the 
center of cell c. 



following will be discussed in terms of the latter. The 
tensor Q is symmetric and is defined as 

Q = iri(a+C'j)e (3) 

where e — (^j^q) is a 7r/2-rotation matrix and 



(4) 



Here i,j are Cartesian components, v denotes a vertex 
of the cellular solid, the definitions of the vectors Tcv and 
Rev are given in fig. 1, and the sum runs over all the 
cells c that surround vertex v. The only contribution 
to Q comes from the parallel parts of Rev and rev and 
therefore Q is a measure of the net rotation of the triangle 
around vertex v relative to its immediate environment. It 
is also the net deviation of the quadrilaterals formed by 
the self-dual r-R pairs (see fig. 2) from rhombi in a given 
direction. The geometric interpretation of Q can be best 
seen in terms of P = iQe^^ . Using vector algebra, P can 
be rewritten in terms of the vectors Scv and tcv of fig. 2 
as 



(5) 



The antisymmetric part of each term in Cv can be writ- 
ten as Acv^ where Acv is the area of the quadrilateral 
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FIG. 2. The quadrilateral of area Acv whose diagonals are 
fcv and Rev ■ The vectors Scv and tcv appear in the symmetric 
tensor P. Note that tcv = —Scv' and therefore a sum over P 
inside cell c vanishes identically irrespective of the shape of 
the cell. 



whose diagonals are Rev and Tcv The area associated 
with vertex v is then Ay = Acv and the total area of 
the system is exactly Asys — X^t, ^v- 

Although the first-principles derivation of Q encourag- 
ingly supports cq. (2) this equation can only be useful if 
Q can bo coarse-grained to yield a macroscopically valid 
equation. This, however, proved to be non-trivial [1] se- 
riously impeding the application of the new theory. This 
problem is resolved here. 

Coarse-graining of constitutive properties has been 
studied in many contexts. The conventional approach 
is to integrate over small-scale degrees of freedom, leav- 
ing the constitutive properties dependent on large-scale 
coordinates [6] . The procedure is to choose a lengthscale 
L above which the d-dimensional medium appears statis- 
tically homogeneous. The entire system is then regarded 
as made of basic units of volume L'' and the constitutive 
property is averaged within a unit, giving a mean charac- 
teristic value. This value is then used in the constitutive 
equation, which is subsequently regarded as applicable to 
scales larger than L. This straightforward approach fails 
for the constitutive eq. (2) because the mean value of 
the tensor Q vanishes, giving a trivial identity for macro- 
scopic scales. To see this consider a part of the system 
r of area A^ and boundary dV, which consists of 
vertices. The average of the tensor P in F is 



(6) 



Prom fig. 2 we see that if v and v' are neighbor vertices 
then tcv = —Scv' in the cells c that they both border. This 
means that the terms in sum (6) cancel out in pairs for 
cells that are fully enclosed in T and only terms along the 
boundary dT contribute. It follows that Pp diminishes 
with size as 1/ v^iVr- Therefore a simple area average of 



Q vanishes on large scales irrespective of the microstruc- 
tural characteristics. To use eq. (2) for macroscopic sys- 
tems and keep Q as the correct geometrical decsriptor 
we need to develop a coarse-graining method that goes 
beyond straightforward area averaging. 

The key to the resolution of this problem comes from 
the observation [1] that in systems possessing a staggered 
order (stagger-ordered systems, SOS) eq. (2) can be 
rewritten in terms of only half the degrees of freedom. 
A staggered order arises when it is possible to label all 
the vertices -|- and -, such that each vertex has neigh- 
bors of only the opposite sign. This requires that all 
the cells in the network have an even number of edges. 
Such networks can be partitioned into -|-/- vertex pairs, 
each of area Apair = + A~ . The local force mo- 
ment on the triangle around v is Sv = ^^v' x fw', 
where v' are the neighbors of v, fw' the forces that tri- 
angles v' exert on v, and fw' are the position vectors of 
the contacts between v and v' . By definition, the local 
stress at vertex v is av = Sv/Av Let us define the mean 
stress and the stress difference of a pair, respectively, as 
CTm = {S~^ + S~)/Apair and ad = {S^ - S~)/Apair- o'm is 
the stress field, coarse-grained over a scale of two vertices 
and (Jd is the local fluctuation whose average must vanish 
on larger scales. A manipulation of eqs. (2) for both the 
+ and - vertices gives 

(g+ + Q-) : <t™ + (g+ - 0") : <Trf = 
(0+ - Q-) : am + (0+ + Q-) ■■ cTd = 

In an area average of these equations the term linear in 
(Q^ + Q^) vanishes, as discussed above, and ad van- 
ishes from its definition. This renders the first equation 
a trivial identity. It also means that the second term 
in the second equation consists of a correlation of two 
small quantities and is therefore negligible relative to the 
first term. Thus only the first term of the second equa- 
tion survives the averaging. Furthermore, the vanishing 
of the average of + Q~ means that the average of 
— Q~ is equal to that of 2(5"'". Thus, on a lengthscale 
of several cells 
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This relation has the same form as the original consti- 
tutive equation and is therefore a coarse-grained version 
of it. Its advantage is that, unlike Q, an area average of 
(5+ does not vanish identically [7] and so the problem of 
coarse-graining is resolvable in SOS. 

Unfortunately, most foams are not SOS and cannot 
be conveniently partitioned into -\-/- vertex pairs. The 
staggered order breaks around cells with odd numbers 
of edges (odd-edged cells, OECs) which cannot perfectly 
accommodate alternating signs around them. At these 
locations the system is 'frustrated' much like a system 
of antiferromagnctically coupled Ising spins. To resolve 
the coarse-graining problem for general systems we then 
need to address the issue of frustration. 
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FIG. 3. The frustration around an OEC (shaded) sur- 
rounded by even-edged ceUs. Labehng the vertices around it 
alternately by -I- and - leaves two neighboring 'frustrated' ver- 
tices with the same sign {v and «'). Starting from a neighbor 
of V and labeling the next shell (thick lines), leaves another 
frustrated pair lying adjacent to the first one. Continuing 
this procedure results in a line of such frustrated vertex pairs 
emanating from the OEC, the first three of which are shown. 
This line can be capped only by encountering another OEC 
along its path. 

Let us consider first one OEC embedded in a region 
of cells with only of even numbers of edgeds (e.g. fig. 
3). Starting from an arbitrary vertex, v, we label the 
vertices alternatingly + and - around the OEC in the 
clockwise direction. The last vertex, v' , has the same 
sign as the first vertex v. Apply now the same procedure 
to the first shell of cells surrounding the OEC, starting 
from the neighbor of v and going clockwise, as shown in 
fig. 3. Since there are only even-edged cells around the 
OEC then this shell must also contain a frustrated pair 
of vertices and, by construction, this pair lies adjacent 
to v-v'. Repeating this process shell by shell outwards 
results in a line of frustrated pairs of vertices emanating 
from the OEC. Three pairs along this line are shown in 
fig. 3 highlighted by ellipses. This line of frustrated 
pairs extends to the boundary of the system unless it 
encounters another OEC along its path, in which case it 
terminates. 

This suggests a procedure to isolate the frustrated 
pairs in general systems and remove their effect: 1. Iden- 
tify all the OECs in the network; 2. partition the OECs 
into nearest pairs; 3. following the above procedure, draw 
a line of frustrated vertex pairs between each pair of 
OECs avoiding intersection between lines. Except along 
these lines, the system is ordered, i.e. it can be com- 
pletely partitioned into pairs of +/- vertices. The choice 
of nearest pairs on the second step minimizes the total 
number of frustrated pairs. 

The frustration is removed by renormalization. Each 
frustrated pair of vertices along the frustration lines (say 
V and v' in fig. 3) is regarded as one supervertex. These 



FIG. 4. Fusion of a frustrated pair of vertices renormalizes 
away the contact point between the triangles. The area of 
the renormalized vertex (shaded) is Af = Ai + A2, which is 
the coefficient of the antisymmetric part of the renormalized 
geometric tensor (7/ = C'l + C2- The common vectors s and t 
cancel and therefore the symmteric part of the renormalized 
grain is independent of the vectors that connect the contact 
point to the centers of the cells that surround it. Therefore 
the stress field, that is coupled to the microstructure via the 
symmetric part of C is unaffected by the fusion. 

vertices (the equivalent of grains) are pressed against 
each other anyway and therefore this does not alter the 
original external forces on the pair. The renormalized 
pair geometric tensor is defined as Cpair ^ Cy + Cv> and 
hence Ppair — Pv + Pv' ■ Since this tensor determines the 
stress through eq. (2) it remains to show that this opera- 
tion leaves the stress field outside the supervertex intact. 
Using eq. (5) and inspecting fig. 2 note that the vectors 
Scv' and tcv, which connect the v-v' contact point to the 
center of cell c, cancel out in this sum. So do their coun- 
terparts in the opposite cell and so Ppair is independent 
of this contact. This also means that the stress tensor be- 
comes independent of fvy> . The remaining s and t vectors 
surround the supervertex as if it were a regular vertex and 
the geometric tensor outside the supervertex remains un- 
affected. Therefore the renormalization operation ineed 
leaves the stress field outside this pair unchanged, as re- 
quired. Another convenient feature of the procedure is 
that the area of the supervertex is Apair = Ay + A^i, 
exactly corresponding to the sum of the adjacent areas 
of the two original grains (see fig. 4) and so the volume 
of the system is preserved. In granular packings this pro- 
cedure is the same as actual fusion of the pair of grains. 
Once all the frustrated pairs are renormalized the system 
becomes fully SOS and the above coarse-graining can be 
used. This procedure then makes it possible to upscale 
the constitutive eq. (2) in general systems. 

To illustrate the method, consider the structure shown 
in fig. 5. It is a honeycomb lattice in which one cell 
has been made into an octagon at the expense of two 
neighbors that are reduced to pentagons. Except around 
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FIG. 5. A honeycomb cellular structure containing a defect 
which consists of an octagonal cell straddled by two pentago- 
nal cells. The coarse graining is over all the positive vertices 
but only those near the defect with irregular areas (shown 
shaded) have a finite P. 



the pentagons the structure is SOS with P vanishing 
due to the hexagonal symmetry. The structural defect 
gives rise to a fluctuation in the constitutive proper- 
ties, which we wish to determine. Following the pro- 
cedure outlined above, we renormalize the two pentag- 
onal cells by fusing the two pairs enclosed by ellipses. 
Summing over the Q-tensors of the positive vertices gives 

Qdefect = ( ^4v^ I ) and, after some algebra, the con- 

5 2 

stitutive equation reads 

16 

(Til + 5fT22 + -^cri2 = (8) 
v3 

Combined with eqs. (1), this relation gives the stress on 
the scale of the entire defect. 

To conclude, this communication has addressed the ap- 
plication of the microscopic theory for stress transmis- 
sion in trivalent solid foams and marginally rigid gran- 
ular assemblies to macroscopic lengthscales. Under con- 
ventional coarse-graining the equation that couples the 
stress field to the microstructure becomes uselss an up- 
scaling method has been proposed to overcome this prob- 
lem. The method consists of fusion of frustrated ver- 
tices / grains, leading to a renormalized system that is 
fully stagger-ordered and whose structure is character- 
ized by a suitably renormalized geometric tensor. The 
renormalized system can be partitioned into two subsys- 
tems, much like antiferromagnetic sublattices, and this 
enables coarse-graining of the constitutive equation. The 
method makes it possible to extend the microscopic the- 
ory to macroscopic systems, paving the way for calcula- 
tions that are useful for engineering purposes. An ex- 
ample of a localized defect in an otherwise perfect hon- 
eycomb lattice; has bcicn solved explicitly. This solution 
can be used to give the Green function for the stress field 
in generally disordered lattices but this issue will be dis- 



cussed in more detail elsewhere. Since the microscopic 
theory does not resort to stress-strain relations then nei- 
ther does the upscaled macroscopic version. This obvi- 
ates the problem of the structure-property relationship as 
we know it in these systems. Specifically, there is no need 
to take the conventional approach where a relationship 
is sought between the local structure and elastic con- 
stants, which are subsequently coarse-grained. Rather, 
the problem reduces to the direct coarse-graining of the 
structural tensorial descriptor Q. The new method also 
makes it straightforward to predict stress concentration 
in cellular materials and relate the occurrence of 'hot 
spots' to local microstructural defects. This is expected 
to improve the prediction of failure in advanced mate- 
rials as well as to set new goals for material processing 
methods to manipulate microstructural characteristics in 
order to avoid such vulnerability. 
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